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We consider the equiprobable distribution of spanning trees on the square lattice. All 
bonds of each tree can be oriented uniquely with respect to an arbitrary chosen site called 
\Q | the root. The problem of predecessors is finding the probability that a path along the oriented 

bonds passes sequentially fixed sites i and j. The conformal field theory for the Potts model 
predicts the fractal dimension of the path to be 5/4. Using this result, we show that the 
• probability in the predecessors problem for two sites separated by large distance r decreases 

as P(r) ~ r -3 / 4 . If sites i and j are nearest neighbors on the square lattice, the probability 
= 5/16 can be found from the analytical theory developed for the sandpile model. The 
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c*3 . known equivalence between the loop erased random walk (LERW) and the directed path on 
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the spanning tree says that -P(l) is the probability for the LERW started at i to reach the 
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neighboring site j. By analogy with the self-avoiding walk, -P(l) can be called the return 
probability. Extensive Monte-Carlo simulations confirm the theoretical predictions. 

E-mail: vahagn.poghosyan@uclouvain.be, priezzvb@theor.jinr.ru 



■ Keywords: Loop erased random walk, spanning trees, Kirchhoff theorem, Abelian sandpile model. 

in 

IT) ■ I. INTRODUCTION AND MAIN RESULTS 
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In the graph theory, the spanning tree of connected graph G is a connected subgraph of G 
containing all vertices of G and having no cycles. Numerous applications of spanning trees began 
with the seminal Kirchhoff's problem solved in 1847 and then spread out many branches of mathe- 
matics and theoretical physics. In the statistical mechanics, spanning trees are related to the Potts 



model [l| , the dimer model 



2]. the sandpile model [jj and many others. A relation between lattice 



models of statistical mechanics and spanning trees via the Tutte polynomial has been established 



by Fortuin and Kasteleyn 
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The Kirchhoff theorem claims that the number of spanning trees of connected graph G is a 
cofactor of the Laplacian matrix A of graph G. If one deletes any row and any column from A, 
one obtains a matrix A* which gives the number of spanning trees as det A*. The determinantal 
structure allows easy calculation of local characteristics of the spanning trees, for instance, the 
average number of vertices with given number of adjacent bonds. A characterization of non-local 
objects in the spanning tree is not so simple. One of such the objects is the chemical path defined 
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as a path along two or more bonds of the tree. The fractal dimension of a long chemical path 
on the two-dimensional lattice has been calculated by means of a mapping of the spanning tree 
configurations onto the Coulomb gas model. 

A closely related object is the loop erased random walk (LERW) on the two-dimensional lattice 



111 ] which was proven to be equivalent to the directed chemical path of the spanning tree of the 
same lattice In this paper, we consider a problem arising in the theory of LERW and 

equally distributed spanning trees: given two lattice sites i and j, what is the probability that the 
LERW or the directed chemical path passes i and j. If site i is passed first, we say that i is the 
predecessor of j and coin the mentioned problem as the predecessor problem. Surprisingly, the 
problem has no exact solution in the general case. Only two limiting cases are available: (a) If 
sites i and j are separated by large distance r , the asymptotics of P(r) can be found from known 
results on the fractal dimension of the chemical path; (b) If points i and j are nearest neighbors of 



the square lattice, the seeking probability can be found from the theory of sandpiles |9j (see also 

The asymptotic behavior of P(r) for large distance r follows directly from the definition of 
fractal dimension. Indeed, consider a large square lattice C and the set of uniformly distributed 
spanning trees on C We assume that the root is situated at the boundary of C. Consider site i 
in the bulk of the lattice and some circle contour C of radius R with the center in i. Let LT be 
a directed chemical path from i to the root along the oriented bonds of a tree. All points of the 
subset of LT inside C are descendants of i. In accordance with the definition of the fractal dimension 
of the directed path on the spanning tree, the number of the descendants inside C is proportional 
to i? 5 / 4 (see Majumdar [7]). The probability that point i is the predecessor of point j lying at 
distance r from i is the density of descendants p{r). Thus, we have 



pR 

J p(r)rdr ~ i? 5 / 4 (1.1) 

from where we conclude that p(R) ~ R^ 3 ^. 

As it was mentioned before, the problem of predecessors for an arbitrary disposition of two 
lattice points is not solved. In Section 2 we concentrate on a particular problem of probability 
-P(l) when points % and j are nearest neighbors of the square lattice. An essential element of the 
theory of sandpiles is the probability distribution of sites having 0, 1, 2 and 3 predecessors among 
the nearest neighbors. The corresponding probabilities are denoted by Xq, X±, Xi and A3. Having 
explicit expressions for these values, we obtain P(l) as their combination and get an unexpectedly 
simple result P(l) = 5/16. In Section 3 we relate this result to the return probability of the LERW. 
Section 4 contains results of the Monte-Carlo verifications. 
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II. THE PROBLEM OF PREDECESSORS FOR NEAREST NEIGHBORS 

The spanning tree enumeration method, namely, the Kirchhoff theorem is proved to be a pow- 
erful mathematical instrument for the investigation of various combinatorial problems of the theo- 
retical physics. In the last decade, it has been developed and adapted for the calculation of height 
probabilities of the Abelian sandpile model [j], [lfj. The Abelian sandpile model is a stochastic 
dynamical system, which describes the phenomenon of self-organized criticality. During the evolu- 
tion the system falls into a subset of all possible states, called the subset of recurrent states. The 
problem is to calculate analytically various observable values in the recurrent state, such as height 
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probabilities Pi, i = 1, 2, 3,4 at a fixed site and height correlations between distinct fixed sites 

It was shown that the calculation of height probabilities in the Abelian sandpile can be reduced 
to the calculation of Xq, X\, X2 and A3 in the spanning tree model. The exact relation between 
these quantities is given by 

p i = X ; P2 = Pl + ir ; P 3 = ^2 + ^5 P4 = P 3 + X 3 . (2.1) 

Majumdar and Dhar [5] have found in 1991 the probability of height 1, constructing the cor- 
responding defect lattice for the situation when a site io has no predecessors and calculating the 
determinant of the defect matrix A*. A technique for computing the numbers X±, X2, A3 has been 
devised in j^. The results are (see also for details) 

8(tt -2) 3 48 15 3 1 48 6 3 16 1 3 

Xo — 3 — ; Ai — - + ^ 2~n~ ' ^-2 — 7 3 + ~ 2+ - ! -^3 — - 3 H — 2 — o - ' v 2 - 2 ) 

7T- 3 4 7T J 7T Z ZTT 4 7T- 3 7T Z 7T TT 7T Z ZTT 

which give 

2(vr- 2) 1 12 3 1 3 12 1 3 4 1 1 

P\ = 5 ; "2 — 7 + ^ o ~~ o~~ ' " 3 ~~ a 3 + _ ' " 4 ~~ o + ~~ 3 + ~~2 ~~ o~~ ■ ( 2 ' 3 ) 

TT 6 4 TT TT ZTT 8 TT A TT 8 TT TT Z ZTT 

Now consider the problem of predecessors for nearest neighboring sites. First we fix a site io 
in the bulk of the square lattice. Denote its right nearest neighboring site by jo- Next consider 
four various cases, when io has exactly k nearest neighboring predecessors (k = 0,1,2,3) (see Fig 
[T]). For k = the site jo trivially is not a predecessor of io- For k = 1, we have 1 of 4 equivalent 
situations when jo is the predecessor of io- For k = 3, we have 3 of 4 equivalent situations when 
jo is the predecessor. The crucial case is the k = 2. Here we have 6 situations, but not all of them 
are equivalent. On the other hand, we can select two groups of 4 and 2 elements (the first four 
and the last two in the third line of Fig. 1) so that elements in each group are equivalent. We are 
looking for the situations where jo is a predecessor of io- There are 2 encircled elements from the 
first group and one from the second group, which correspond to the desired situations. Thus, if we 
take the linear combination of X\, A2 and A3 with coefficients 1/4, 1/2 and 3/4 correspondingly, 






FIG. 1: All possible situations with a fixed vertex io (the central vertex in diagrams) to have various nearest 
neighbouring predecessors on square lattice. By black color we denote vertices, which are predecessors of 
io. White means that the corresponding vertex is not a predecessor of io- On the k + 1st row (k — 0, 1, 2, 3) 
situations with k predecessors are presented. By the dashed line are circled the configurations for which the 
right neighboring vertex is a predecessor of io- 

we get the desired probability -P(l) that jo is the predecessor of io: 

P(l) = -X 1 + -X 2 + -X 3 = — . (2.4) 
v ; 4 2 1 4 16 K ' 

III. RETURN PROBABILITY FOR THE LOOP ERASED RANDOM WALK 

Consider a finite square lattice C with vertex set V and edge set E. Given V = [uo,ui,U2, ■ ■ ■ , u n ] 
a path in £, its loop-erasure LE('P) = [70,71,72, ••• ,7m] is defined by chronologically removing 
loops from V . Formally, this definition is as follows. We first set 70 = uq. Assuming 70, ... ,7fc 
have been defined, we let Sk = 1 + max{j : Uj = 7^}. If si- = n + 1, we stop and LE("P) = 
[70,71,72, •• • ,7m] with m = k. Otherwise, we let 7fc+i = u Sk . Note that the order in which we 
remove loops does matter, and it follows from the definition, we remove loops as they are created, 
following the path. A walk, obtained after applying the loop-erasure to a simple random walk 
path is called Loop-Erased Random Walk (LERW) . Since the infinite simple random walk on finite 
connected undirected graphs is recurrent, the infinite LERW is not defined. On the other hand, 
we can fix a subset W C V of vertices and define LERW from a fixed vertex uq to W . To do that 
we take a path of a simple random walk started at uq and stopped upon hitting W, after that we 
apply loop-erasure. 

Wilson established an algorithm to generate uniform spanning trees, which uses LERW 131 ] . It 



turns out to be extremely useful not only as a simulation tool, but also for theoretical analysis. It 
runs as follows. Pick an arbitrary ordering V = {vq, v\, . . . , vn} for the vertices in C. Let Sq = vq. 
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Inductively, for i = 1,2, . . . , N define a graph Si to be the union of Si_i and a (conditionally 
independent) LERW path from Vi to £j_i. Note, if Vi G Sj_i, then = £j-i. Then, regardless of 
the chosen order of the vertices, SV is a UST on £ with root vq. If we take as an initial condition 
So = W, with some W C V, then the generated structures will be spanning forests with set of 
roots W. The spanning forest with fixed set of roots can be considered as a spanning tree, if we 
add an auxiliary vertex and join it to all the roots. 

Now consider the Wilson algorithm on C with the set of boundary vertices dC and take So = dC. 
When the size of the lattice tends to infinity, the boundary effects will vanish, so we can neglect the 
details of the boundary. So we will not distinguish between spanning forests and spanning trees. It 
follows from the Wilson algorithm for a fixed site io, that if we start a LERW from iq upon hitting 
the boundary dC, we will generate a path £ of a spanning tree from io to the boundary (see also 
Q], sjj). All vertices on the path t form the set of descendants of io- So, if a fixed vertex jo belongs 
£, iq is a predecessor of jo- 

Consequently, the probability P(jo — io) that io is a predecessor of jo in a randomly taken (from 
the uniform distribution) spanning tree on the large square lattice equals to the probability that the 
LERW started from io passes jo- In the particular case |jo — io| = 1, the probability -P(l) = 5/16 
calculated in the previous section is the return probability for the LERW. 



IV. MONTE-CARLO SIMULATIONS 



Consider finite 2N + 1 x 2N + 1 square lattice £j\r. Denote its central vertex by io an d assume 
that it is an origin of the coordinate system. We deliberately took an odd-odd lattice to provide the 
symmetry which guarantee the most efficient vanishing of boundary effects for large lattices. After 
generating a large amount of LERWs starting from io, we get an approximation of the function 
P/v(jo — io) = PnUo)- Given fixed jo, we can extrapolate PnUo), tending N to infinity and get 
asymptotical function P(Jo)- Assume that the Euclidean distance between the origin io and jo 
is r and coordinates of jo are (reos</?, rsimp). The Monte-Carlo simulations and Coulomb gas 
arguments show that the asymptotical behaviour of the function P(jo) for large r (r > 1) does not 
depend on (p. So, for r>lwe have P (jo) — P{r). Fig. [2] shows the behaviour of P(r) for various 
jo on the log-log scale, obtained from Monte-Carlo simulations. From this result we conclude that 
P(r) decreases as a power function: 

P(r) ~ ^, (4.1) 

with a ~ 0.751 and C ~ 0.305. During the simulations, we took sizes up to N = 100 and number 
of simulations 10 8 . The obtained results are in agreement with a = 3/4 which follows from the 
Coulomb gas arguments mentioned above. The effective Monte-Carlo algorithm allows evaluation 
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FIG. 2: The results of Monte-Carlo simulations for the probability P{r). 



of probabilities P(jo — io) for arbitrary finite jo, Iq. At the same time, the analytical calculation of 
P(r) for r > 1 remains a difficult unsolved problem. 
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